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t^- ' Abstract 

O ! 

For the case of small matter effects: V <C Am /2E, where V is the matter potential, 
we develop the perturbation theory using e = 2VE/Am 2 as the expansion parameter. 
We derive simple and physically transparent formulas for the oscillation probabilities in 
the lowest order in e which are valid for arbitrary density profile. The formulas can 



be applied for propagation of the solar and supernova neutrinos in matter of the Earth, 
substantially simplifying numerical calculations. Using these formulas we study sensitivity 
of the oscillation effects to structures of the density profile situated at different distances 
from the detector d. We show that for the mass-to-flavor state transitions, e.g., v<i — > u e , 
the sensitivity is suppressed for remote structures: d > l v E/ AE, where l v is the oscillation 
\£> ■ length and AE/E is the energy resolution of detector. 

o 
^t- . 

§ ! 1 Introduction 

o ; 

For the LMA oscillation parameters the oscillations of solar and supernova (low energy) neu- 

I 



trinos inside the Earth occur in the weak matter effect regime. That is, the matter potential, 



Ph! V, is much smaller than the kinetic energy of the neutrino system: 



where 

C3 



V(x) = V2G F N e (x), (2) 

Gf is the Fermi constant, N e (x) is the number density of the electrons, and Am 2 = m 2 . — m 2 
is the mass squared difference. 

In this case one can introduce a small parameter 

, x 2EV(x) 

and develop the perturbation theory in e(x). 

Neutrino oscillations in the weak matter effect regime have been discussed extensively before 
[1] - [20], in particular, for the solar and supernovae neutrinos propagating in the matter of the 
Earth. The previous work has been done mainly in the approximation of density profile which 
consists of several layers with constant density. 

Here we derive general formula which is valid for arbitrary density profile provided that the 
condition (1) is satisfied. The formula is the generalization of our result in [20] obtained for 
thin layers of matter. 
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2 e- perturbation theory 



We will consider two active neutrino mixing 

v f = U(9)is mass , (4) 

where Vf = (v e , u a ) T ', v mass = (z/i, u 2 ) T are the flavor and mass states correspondingly, and 

_ / cos 9 sin9 \ , . 

~ y — sin 9 cos 9 J 

is the mixing matrix in vacuum. In general v a is a combination of z/ M and v T . 

We will first find the evolution matrix for the mass eigenstates and then make projection 
onto the flavor states. Evolution of the mass states is given by the equation 

.dv mass , , 
i dx = H{x)v mass (6) 



with the Hamiltonian 



The Hamiltonian can be rewritten as 



where 

A m (x) = ^^(cos29-e(x)y + sm 2 29 (9) 

is the difference of energies of the neutrino eigenstates in matter u m = (u lm , v 2m ) T ■, and U' = 
U(9') is the instantaneous mixing matrix of the mass states in matter with the angle 9'{x) 
determined by 

. dx) sin 29 , , 

sin2^(x) = . v ; (10) 

v /(cos2^-e(a;)) 2 + sin 2 2^ 

Using (9) the expression for the mixing angle of the mass states in matter, 9'(x), can be rewritten 

as 

. n „. Vsin2^ . n „ 

sin 20' = — = esin20 m . (11) 

In the last equality of Eq. (11) 9 m is the mixing angle of the flavor states in matter. Taking 
into account relations Vj = U{9 m )u m , and u mass = U{9')u mi we find that Vj = U m {U')^v mass , 
and consequently, U = U m U'. Therefore 

9' = 9 m - 9. (12) 
The angle 9' is small: sin 29' < e, as it should be in the weak matter effect regime. 

The formal solution of the equation (6), that is, the evolution matrix from the initial point 
x to the final point Xf, can be written as 

S(x ^ x,)=Te- i % Hlx)dx , (13) 
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where T means the chronological ordering. Let us divide a trajectory of neutrinos into n equal 
parts (layers) of the size, Ax, so that n — (xf — Xq)/ Ax, and assume constant density inside 
each layer. Then for the evolution matrix we obtain 

S(xo — > Xf) = e~ iH ( Xn ^ Ax ■ e ~ iH ( x ^-i) Ax . . . e ~ iH ( x j) Ax . . . e ~ iH ( x i) Ax 
(x n = Xf). According to (8) it can be written as 

S(x - x f ) = U' n D n U* ■ CdAi-ittf-i ■ ■ ■ UjDjU? ■ ■ ■ U'^U? , (15) 

where 

D i=[l ) > W = ^{V 3 )Ax (16) 

is the evolution matrix of the matter eigenstates in j-th layer, 0™ is the relative phase between 
the matter eigenstates acquired in the layer j, and Vj is the value of potential in the j-th layer. 
The mixing matrix C/j is given by 

,_ / cos^- sin^. \ 
Ui \ -sine'j coser J ' 1 > 

where 6j is the mixing angle of the mass states in the layer j. 

Structure of the expression (15) is rather transparent: it is the product of the blocks 
{U'jDjU'j) for all the layers. In each layer we first project the mass states onto the matter 
eigenstates u m , then evolve the eigenstates and then project back to the mass states. 

Using (16) and (17) we find expression for the j-th block: 

U'DjU? = Dj + G 3 , (18) 



where 



G 3 = (e^ - l) 



1 • n m ( 1 \ . 2 „/ M 
2 Sm2 M 1 + sm M -1 



(19) 



Notice that Gj = 0(e), furthermore the second term in (19), being proportional to sin 2 9', is of 
the order e 2 . 

Inserting expression (18) into (15) we find 

S(x -> x n ) = (D n + G n )(D n _ 1 + G n _0 •••(£>! + d), (20) 

and in the form of series in powers of G — 0(e) it can be written as 

S { x -> x f ) = D n D n _ x ■■■D 1 

n 

+ J2 D n---D i+1 G j D j _ 1 ---D 1 

3=1 
n j-1 

+ E E D n ■ ■ ■ Dj+iGjDj-i ■ ■ ■ D k+1 G k D k -! ...£>! + ..• (21) 

j=l k=l 
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The products of Dj which appear in this formula equal 



n j=fc ...n^- = exp i E <t>7 =exp i J2 &m(V 3 )Ax . (22) 

\ j=k...n } \ j=k...n ) 

In the limit n — > oo and Arc — > the sums are substituted by the integrals: J2 Ax — > / dx. 
So that 

n i =*... w D J --exp(<^^ eB ), (23) 



where 



Furthermore, 



C fc ->*„ = / dxA m (x). (24) 
Gj -> iT(x)da:, (25) 



where according to (19) 



T = 1 sin 2(9 ( J J j + A m (x) sin 2 6K(x) ( J ° J . (26) 

Here we have taken into account the relation (11). 

Substituting the sums by the integrals in Eq. (21) we find the S-matrix in terms of T: 





e^o 



-II dX L dy {l e^) T{x) {l e4J T ^(j e 4J+"--( 27 ) 



Essentially we perform an expansion of the .S-matrix in powers of T. 

Let us take the zero and the first order terms in T in (27) which correspond to the zero and 
first order terms in e. We find 

, / 1 \ .sin 2^ rxt ( e^o-A 

S(*o - *,) = [ Q + .— X o cfaV(x) ^ eW _ B/ Q j . (28) 

An interesting feature of the formula (28) which leads to important consequences (see sec. 
4) is that in the second term the 12-element is determined by the phase acquired from the 
initial point Xq to a given point x, whereas the 12-element depends on the phase from a given 
point x to the final point of evolution Xf. 

The key feature of the method is that we calculate the adiabatic phase exactly, whereas the 
mixing and the amplitude of oscillations are found approximately using the perturbation in e. 

Notice that the second term in T (26) (of the order e 2 ) is proportional to the diagonal 
matrix. The 11-element of this term does not contain oscillatory factor and therefore, being 
inserted in (27), produces the contribution to the S- matrix which is proportional to the length 
of the trajectory. It can be shown that this large contribution is canceled precisely by the 
contribution from the term of the second order in T (second line in Eq. (27)). 
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3 Probabilities of neutrino conversion 

Using the evolution matrix in the mass states basis (28) we can calculate the amplitudes and 
probabilities of different transitions. The evolution matrix in the flavor basis, Sf, can be 
obtained immediately: 

S f = USUI (29) 
The evolution matrix from the mass states to the flavor states, Sf m , is then 

S fm = US, (30) 

and U is the vacuum mixing matrix (5). 

Let us consider the most important examples. The amplitude of the mass-to-flavor transi- 
tion, z/j — > v a , on the way from x to Xf is given by 

A-v^Va = U a jSji. (31) 

Inserting the matrix (28) in this expression we find the amplitude of the v\ — > v e oscillations 
in the first order in V: 

1 f' x f 'A*m 

A Ul ^ Ue = C os9 + - sin 29 sin dx V{x)e^*f . (32) 

2 Jxo 

Then the probability of the v\ — > v e oscillations equals 

P ui ^ Ue = cos 2 9-\ sin 2 29 P dx V{x) sin 0™ . (33) 
2 Jx 

The probability of v 2 — > v e transition relevant for the solar neutrino oscillations in the Earth 
can be obtained immediately from the unitarity condition: 

P U2 ^ Ue = sin 2 9 + I sin 2 29 P dx V(x) sin 0™. (34) 
2 Jx 

Then the regeneration parameter defined as f reg = P U2 ->u e — sin 2 9 (see, e.g., [10]) equals 

f reg = ^sin 2 2fl P dx \ '(./•) >in o"' (35) 
Using expression for the phase 4> x n L^ Xf (24) we obtain explicitly 



f reg = ± sin 2 2£ jT cfe sin jT */^^ (cos 2^ - + sin 2 20. (36) 

Let us underline that these expressions are valid for arbitrary density profile provided that 
the condition (1) is fulfilled. 

The terms which contain V are implicitly of the order e. Indeed, noting that the phase 
4> x ^x f is proportional to Am 2 /2E and changing the integration variable x — > xAm 2 /2E, one 
obtains the factor e(x). 
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It is straightforward to show that the obtained formulas reproduce the well know results for 
particular density distributions. Thus, for one layer with constant potential the Eq. (35) gives 
immediately: 

f reg = e sin 2 29 sin 2 ^ , (37) 

I'm 

where l m = 2n/A m is the oscillation length in matter. For the profile with n symmetric shells 
we obtain performing explicit integration in (35) 

2E n ~ l 

freg = ~T 2 ^ 26 sil1 $ E AV i Sil1 $ * > ( 38 ) 

Am i=Q 

where $o is the total phase acquired along the trajectory and $; is the phase acquired inside 
external radius of a given shell, AVi is the jump of potential at border of the i-th shell [19]. 
In general, the formulas can be obtained from the adiabatic perturbation theory, performing 
integration by parts in Eq. (18) of [19]. In the lowest order in e the adiabatic perturbation 
theory and the e-perturbation theory coincide. 



4 Sensitivity of oscillations to structures of the density 
profile 

The formulas obtained in section 3 can be used to estimate a sensitivity of the oscillation effects 
to various structures of density profile. 

Suppose that there is some structure in the density profile (or V) in the point x. Then 
according to (33), (34) or (35) this feature will be integrated with sm<f) x _, Xf - the periodic 
function whose phase is acquired from the structure position to the detector: d = Xf — x. The 
larger the distance and therefore the phase, the stronger averaging effect is expected. So that 
the effect of remote structures of the profile on the mass-to-flavor oscillation probabilities is 
suppressed. 

To quantify this sensitivity let us introduce the energy resolution function f(E', E) and 
perform averaging of the probability folded with /: 



p 

1 V\— >u e 



JdE'f(E',E)P Ul ^ e . (39) 



It is convenient to parameterize the effect of integration introducing the attenuation factor F(d) 
in the probability as 



| rXf 

P vi ^v a = cos 2 6> sin 2 20 I dxV(x)F(x f - x) sm<f) x ^ Xf , (40) 

2 JXQ 



so that in the absence of averaging F — 1. 

Let us take for simplicity the box like resolution function f(E', E) for which 

_ 1 r E +^f 

P — = AEL^ dEP —- (41) 
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Then assuming that AE <C E and making use of the approximation A m ~ Am 2 [l — e cos 26] /2E 
we find 

F(d) = ^LsinQ(d), Q(d) = !^, (42) 

where Z„ ~ / m is the oscillation length (see Fig. 1). 

As follows from the Fig. 1, the factor F(d) is decreasing function of the distance from 
detector. The decrease of F means that contributions from the large distances to the integral 
(40) are suppressed. According to Fig. 1 which corresponds to AE/E = 0.2, the effect of 
structures at distances above 1500 km are attenuated by F at least by factor 5 in comparison 
with those situated near the detector. Correspondingly, the sensitivity to remote structures is 
much weaker. 

The larger AE, the smaller the width of the first peak of F(d). As follows from (42), the 
width and therefore the region of unsuppressed contributions due averaging are given by 

d<l v ^=. (43) 
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This attenuation effect allows us to explain some features of the zenith angle distribution 
of the regeneration factor. No significant enhancement of the regeneration has been found 
in numerical calculations for the trajectories which cross the core of the Earth in comparison 
with the trajectories crossing the mantle only, in spite of the fact that density of the core, and 
therefore e, are 2-3 times larger (see e.g. [10], [11], [13]). The explanation is straightforward: 
the border of the core is at about D = 3000 from the surface of the Earth. Therefore, according 
to Fig. 1 the effect of the structures from such a distance is suppressed by the factor 3-5 
which compensate the increase of e. On the other hand, the effect of even small structures near 
the surface of the Earth can be substantial. 

The conclusions we made so far are valid for the oscillations of the mass states to the flavor 
states. The situation is different for the inverse - flavor-to-mass states transitions. Indeed, for 
v e — > v\ the probability P Ve -+ Vl has similar expression as in (33) but with phase acquired on 
the way from the initial point xq to the position of the structure x: <p XQ ^ x . Therefore with 
the flavor-to-mass transition probability one would probe structures at the opposite (to the 
detector) side of the Earth. This general consideration is in agreement with our results for thin 
layers of matter [20] . 

For the flavor to flavor transition we obtain the probability 

6 m I r xf 

P v ^ Ve = 1 - sin 2 29 sin 2 J^L - - sin 2 29 cos 29 / dx V(x) [sin (j>™ + sin <f% ] (44) 

2 2 Jx 1 

which is sensitive to structures near the surface from the both sides of the Earth. 

The detailed analysis of the averaging effects and other applications of the obtained results 
will be given elsewhere [21]. 

5 Conclusion 

We have derived expressions for the oscillation probabilities in matter with arbitrary density 
profile in the weak matter effect regime: V <C Am 2 /2E. An accuracy of these expressions (33, 
34, 35) is determined by the parameter e. 

The results can be applied for the solar and supernova (low energy) neutrinos crossing the 
Earth. They substantially simplify numerical calculations of the oscillation effects. 

The obtained formulas reproduce the known probabilities for particular density distributions 
(one layer with constant density, several layers, etc.). 

The formulas have very simple structure and can be used efficiently for analysis of various 
effects. In particular, the sensitivity of the oscillation probabilities to various structures of 
the density profile can be easily evaluated. We have found that the mass-to-flavor transition 
probabilities are sensitive to structures situated close enough to detector. Effect of the remote 
structures is attenuated. The distance which can be viewed by a detector is determined by the 
oscillation length divided by the energy resolution of the detector: l u E/AE. These results can 
be used in future in the oscillation tomography of the Earth. 
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